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We study scalar, pseudoscalar, vector, and axial-vector mesons with non-strange and strange 
£Nj " quantum numbers in the framework of a linear sigma model with global chiral U(Nf)z x U(Nf)n 

7-H ' symmetry. We perform a global fit of meson masses, decay widths, as well as decay amplitudes. The 

f*"^ . quality of the fit is, for a hadronic model that does not consider isospin-breaking effects, surprisingly 

C^l ■ good. We also investigate the question whether the scalar qq states lie below or above 1 GeV and 

find the scalar states above f GeV to be preferred as qq states. Additionally, we also describe the 
. axial-vector resonances as qq states. 

1 PACS numbers: 12.39.Fe, 12.40.Yx, 14.40.Be, 14.40.Df 

Cl| I. INTRODUCTION 

Understanding the meson mass spectrum in the region below 2 GeV is one of the fundamental problems of QCD. 
While the quark model seems to work very well for many resonances [see, for instance, the summary in Ref. [l|], 
some fundamental questions, such as the constituent quark content of scalar and axial-vector resonances, are still 
t-H 1 unanswered. 
> 

The quark content of the scalar mesons has been a matter of debate for many decades The Particle Data 
Group (PDG) Q suggests the existence of five I{J PC ) = 0(0++) states in the region below 1.8 GeV: / (500), /o(980), 
/o(1370), /o(1500V, and /o(1710). Also the existence of a sixth state, /o(1790), close to, but distinct from, /o(1710) 
has been claimed @. Additionally, there are scalar resonances in the isotriplet and in the isodoublet sector as well: the 
1=1 states ao(980) and eto(1450) are well established; the existence of the 1=1/2 state ^(800) (or k) is confirmed 
by some Q and disputed by other authors Q, whereas the I = 1/2 scalar kaon, Kq (1430), is an established resonance. 
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A description of all mentioned scalar states as qq states is not possible for a simple reason: the number of physical 
resonances is much larger than the number of resonances that can be constructed within a qq picture of mesons. 
. Explicitly, there is only one I(J PC ) = 0(0 ++ ) state that can be constructed from the (non-strange) u and d quarks 
$_i ' (provided they are degenerate). Denoting this state as er/v we obtain on = (uu + dd)/^/2. An additional I(J PC ) = 
0(0 ++ ) state can be constructed, if the strange quark is considered as well: the pure strange state erg = ss. Since o^ 
and as have the same I(J PC ) quantum numbers, we expect the physical spectrum to consist of mixed, rather than 
pure, states of on and erg; however, such mixing will, of course, produce exactly two new states. These states will 
correspond to at most two of the mentioned five (six) /o states, and the natural question is then: which two? 

Similarly, in both the 1=1 and 7=1/2 sectors one can construct only one quark-antiquark resonance. Restricting 
for example to electric charge +1, one has the states clq = du and = su. The state = du can be assigned to 
ao(980) or ao(1450), and Kq + = su to Xq(800) or i^Q(1430). The question is: which is the correct assignment? 

An answer to these questions is inevitably complicated for several reasons. Firstly, as already indicated, states 
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with the same I(J PC ) quantum numbers are expected to mix - this needs to be considered in particular in the 
scalar sector due to the large number of physical resonances. Secondly, basic features (pole mass and decay width) 
of some scalar mesons [e.g. /o(500)] are notoriously difficult to determine experimentally, which makes it non-trivial 
to determine the structure of these states [see, e.g. Ref. Q for an example regarding the /o(500) resonance and Ref. 
0| for an example regarding / (1370)]. 

The question of how to correctly determine the quark structure of the mesons is not only restricted to the scalar 
sector. More recently, also the nature of the axial- vector mesons, most notably ai(1260), but also that of the isoscalar 
states /i (1285), /i(1525), and the kaonic state i^i (1270) [or A"i(1400), see the discussion in Ref. Q and refs. therein] 
have been investigated. Should one interpret the isotriplet resonance ai(1260) as a quark-antiquark state, as the 
quark model suggests (e.g. = du), or is this state a broad pit molecular state? 

Understanding these issues is not only crucial for hadron vacuum spectroscopy but is also important at nonzero 
temperatures and densities, because the correct identification of the chiral partner of the pion and of the p 
meson is necessary for a proper description of the in- medium properties of hadrons Q. In fact, if the ai(1260) is 
(predominantly) a quark-antiquark state, it is the chiral partner of the p meson, with which it becomes degenerate 
at large T and p The chiral partner of the pion is the scalar-isoscalar state er/v = (uu + dd) /y/2 which has been 
identified with the lightest scalar state /o(500) in many works. However, as various theoretical results for the vacuum 
have shown, and as we shall present in detail in this work, this assignment is not correct: it turns out that /o(1370) 
emerges as the chiral partner of the pion. 

Last but not least, the question of the origin of hadron masses is important. One aims to understand to which 
extent the quark condensate, (qq), and the gluon condensate, (G^ V G^ V ), generate the hadron masses. To give an 
example, which condensate is predominantly responsible for the mass of the p meson? And related to this, how will 
the mass of the p (and that of other resonances) change in the medium? 

The answers to all these fundamental questions are in principle contained in the QCD Lagrangian. Unfor- 
tunately, QCD cannot be solved by analytic means from first principles in the low-energy domain. For this 
reason, effective theories have been developed which share some of the underlying symmetries of QCD. The QCD 
Lagrangian exhibits, in addition to the local SU(3) C color symmetry and the discrete CPT symmetry, a global 
chiral U(Nf)L x U(Nj)r = U(l)y X U(1)a X SU(N/)v x SU(N/)a symmetry which is broken in several ways: 
spontaneously [due to the chiral condensate (qq) = (qnqL + Qt.Qr ) ^ explicitly (due to non- vanishing quark 

masses), as well as at the quantum level [the U(1)a anomaly (l2j)]. 

In the framework of effective theories the chiral symmetry of QCD can be realized along two lines: linearly [l3| 
and non- linearly In the former case one obtains the so-called linear sigma model: both scalar and pseudoscalar 
degrees of freedom are present and vectors as well as axial-vectors can be included into the model in a straightforward 
manner (l5l4T8| . Note that an important advantage of linear sigma models is the possibility to investigate the state 
of matter at large values of temperature and chemical potential [e.g. in the region of the chiral transition, i.e., where 
chiral symmetry is restored (l9l [20|]. Every linear sigma model contains so-called chiral partners - states that mix 
with each other under axial transformations - that become degenerate when chiral symmetry is restored. In the 
non-linear realization (i.e., the non-linear sigma models), the scalar and the (axial-)vector states are integrated out 
and the pseudoscalar states are the only degrees of the freedom: one obtains the Lagrangian of chiral perturbation 
theory [2l| . If the vector mesons are not integrated out, one is left with chiral perturbation theory with vector 
mesons, see e.g. Ref. [22| and refs. therein. 

In this paper, we present a linear sigma model containing scalar, pseudoscalar, vector, and axial-vector mesons with 
both non-strange and strange quantum numbers. Although this project represents a straightforward implementation 
of the principles of the linear realization of chiral symmetry as already outlined in Ref. [l7| , this is - to our knowledge 
- the first time that all these degrees of freedom are considered within a single linear chiral framework. In view of 
the large number of the fields involved, our model shall be referred to as "extended Linear Sigma Model", or "eLSM". 

In Refs. [l5l.l2"3l|. we have already presented a linear sigma model with vector and axial- vector mesons for two flavors. 
Comparing our model with experimental data for meson vacuum phenomenology led us to conclude that the scalar 
qq states are located in the energy region above 1 GeV in the hadron mass spectrum. Explicitly, we concluded that 
the resonances /o(1370) and ao(1450) are strongly favoured to be scalar quarkonia. The present work is more general: 
we now consider three flavors both in the (pseudo)scalar and (axial-)vector channels in order to ascertain whether the 
conclusion of Ref. [l5[ is still valid once strange mesons are considered. We emphasize that, as discussed in Ref. [l5j . 
all fields entering our model describe pure quark-antiquark states. The reason is that masses and decay widths of our 
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theoretical states scale as N® and N^ 1 , respectively, where N c is the number of colors; thus all decay widths vanish in 
the large- N c limit. Note that the inclusion of strange mesons provides us with a large number of very precise data [in 
general, decisively more precise than in the case of non-strange mesons As a consequence, our model parameters 
are much more constrained than in Ref. (la |. Preliminary results of this work have already been presented in Ref. 

In order to (i) test the performance of a linear sigma model in describing the overall phenomenology of (pseudo)scalar 
and (axial-)vector mesons and (it) ascertain which scalar states are predominantly quarkonia, we perform a global fit 
in which 21 experimentally known quantities (decay constants, masses, and decay widths as well as amplitudes) are 
included. The situation in the scalar- isoscalar sector is extremely uncertain because five (six) resonances, /q(500), 
/ (980), / (1370), / (1500), and / (1710) [and, possibly in the future, / (1790)] are listed in the PDG. In addition, the 
/o(500) and / (1370) decay widths and branching ratios are poorly known and /o(980) suffers from a distortion by the 
A A threshold. Moreover, in the present version of the model the scalar glueball is not included, but it is known that 
the latter state with a bare mass of about 1.5—1.7 GeV as predicted by lattice QCD [25[ can sizably affect the masses 
and the branching ratios of the scalar states above 1 GeV, that is /o(1370), /o(1500), and /o(1710) (26[. In view of 
all these reasons, we do not include the scalar-isoscalar states into the fit, but we test a posteriori their assignment 
and phenomenology. We do, however, include the isotriplet and isodoublet quark-antiquark scalar states and we test 
all four combinations a (1450) / AJ(1430), a (980) / A^(800), a (980) / if g; (1430), and a (1450) / A£(800). Quite 
remarkably, the outcome of the fit is univocal: only the pair ao(1450) / Aq(1430) yields a good fit, while the other 
combinations do not. We thus conclude that the 1=1 and 1=1/2 quark-antiquark scalar resonances lie above l 
GeV. In fact, the quality of our fit is surprisingly good. We describe all experimental quantities with an average error 
of 5%, and most of them even to much better precision. We perceive this to be a remarkable achievement within an 
(in principle, quite simple) effective model for the strong interaction. 

We then study the two isoscalar mesons of the model: a consequence of our fit is their mass in the large- N c limit. 
These masses turn out to be about 1.36 GeV and 1.53 GeV, respectively. Varying large- iV c suppressed parameters 
which cannot be determined by our fit, one can also study their decays: the non-strange quark-antiquark state is well 
described by /o(1370) and the heavier ss state might be /o(1710). Turning to the axial-vector channel, both masses 
and decays are well described by assuming the axial-vector resonances ai(1260), /i(1285), /i(1525), and Ai(1270) 
as (predominantly) quark-antiquark states. Detailed calculations of all formulas used in this paper as well as further 
discussion are presented in Ref. (27| . 

The outline of the paper is as follows. In Sec. [IT] we discuss the U(3)l x U(3)r linear sigma model with vector 
and axial- vector mesons. In Sec. Mil we discuss the global fit and its consequences and in Sec. IIVI we present our 
conclusions. We defer detailed calculations to Appendices IA1 and [Bl Our units are H = c = 1; the metric tensor is 
g^ = diag(+, -,-,-). 



II. THE MODEL 



Lagrangian 



In this section we discuss our model: a linear sigma model with (axial-)vector mesons based on the chiral U(3)l X 
U(3)r symmetry. We retain only operators up to fourth order 1 ; higher-order operators cannot appear when dilatation 
symmetry is imposed, see the discussion in Refs. [2^,[2^|. The Lagrangian of the model reads 



C = Tr[(A> Ai $) t (A> Al <I>)] - m^Tr($ t $) - Ai[Tr($ t $)] 2 - A 2 Tr($ t $) 2 
1 



■Tr{L* + R* ) 



Tr 



A (A 



Tr[A($ + &)] 



+ Cl (det $ - det *t f + i^(Tr{A MI/ [L", V] } + Tr{R^ [R" ,R"]}) 

+ y Tr($t$) Tr(L 2 + Rl) + h 2 Tr[(L M $) 2 + ($A M ) 2 ] + 2h 3 Tr (L^R^). 

+ . g3 [Tr(L M L,Z^L l/ ) + Tr^R^R* R u )} + g A [Tr [L^L V L V ) + Tr (R„R^R V R U )} 
+ g 5 Tr {L^) Tr (R V R V ) + g 6 [Tv(L^) Tr(L u L v ) + Tr{R^) Tr{R„R»)] , 



(1) 



1 With the exception of the anomaly term which is of sixth order. 
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where 



= - igi(L»$> - <PR^) - ieA^[T 3 , $] , 
i^' y = d»L v - ieA^[T 3 ,L u ] - {d v L' i - ieA v [T 3l IJ 1 ]} , 
= d»R v - ieA»[T 3 ,R v ] - {d v R' i - ieA u [T s ,R fi }} , 

The quantities P M , and L M represent the scalar and vector nonets: 



* = jT{S i + iP i )T i = -= 



i=0 



i=0 



Rfi = ^L,i V t ~ A i) T i 



i=0 
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f o 

UJN—P _ JIN —O-i 

VI _V2 

A* - A'? 




(2) 



(3) 



(4) 

— /is 

where the assignment to physical particles is shown as well 2 . Here, Tj (i = 0, . . . , 8) denote the generators of f (3), 
while Si represents the scalar, Pj the pseudoscalar, Vf* the vector, and A^ the axial-vector meson fields, and A^ is 
the electromagnetic field. It should be noted that here and below we use the so-called non strange - strange basis in 
the (0 — 8) sector, defined as 



<pz{SuPi,V?,A!t), 



(5) 



which is more suitable for our calculations. Moreover, H and A are constant external fields defined as 

H = PqPo + #8^8 = 



2 








2 U 

n "os 

U V2 



A = A T + A 8 T 8 



^00 
% 

% 



S N 
5 N 
5 S 



(6) 



(7) 



Throughout this work we assume exact isospin symmetry for u and d quarks, such that the first two diagonal 
elements in Eq. © and Eq. ((7J| are identical. 

Let us now discuss the assignment of fields in Eqs. @ - (0]). If we consider isospin multiplets as single degrees of 
freedom, then there are 16 resonances that can be described by the model: er/v, erg, Sq, Aq (scalar); tjn, rjs, 7?, K 
(pseudoscalar); lu^, Wg, p^ 1 . K ni (vector), and f^ N , f^ s , af", K\ (axial- vector) . All fields in our model represent qq 



states, as discussed in Ref. [I5| . If we assign a state from our model to a physical resonance we, therefore, implicitly 
assume that this resonance is a qq state. This assumption can be tested for a multitude of physical resonances in 
the scalar and axial-vector sectors, as discussed below (in the pseudoscalar and the vector channels, there are no 
ambiguities). 

In the non-strange sector, we assign the fields ir and tjn to the pion and the non-strange part of the r\ and rj' mesons, 
T]n = (uu + dd)/V2. The fields and p^ represent the w(782) and p(770) vector mesons, respectively, and the fields 



2 With the exception of the (0 — 8) sector where particle mixing takes place (see below). 
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/Cjy and a/ 1 represent the /i(1285) and ai(1260) mesons, respectively. In the strange sector, we assign the A fields 
to the kaons; the rjs field is the strange contribution to the physical 77 and rj' fields [775 = ss]; the uj§ , fig, A*, and 
A'i fields correspond to the 0(1020), /i(1420), A*(892), and #1(1270) [or Ai(1400)] mesons, respectively. 

Unfortunately, the assignment of the scalar fields is substantially less clear. Experimental data suggest existence 
of five (six) scalar-isoscalar states below 1.8 GeV: /o(500), /o(980), / (1370), / (1500), and / (1710) [as well as 
/o(1790)]. Note that the existence of these five states is acknowledged by the Particle Data Group - PDG [l|. [The 
sixth state, /o(1790), will not be of importance for the rest of our work because its predominant ittt decay renders it 
a putative radial excitation of /o(1370), and our model describes ground-state quarkonia only] 

Our model contains the pure non-strange isoscalar apj and the pure strange isoscalar ag. We will demonstrate 
below that our model yields mixing of <jn and ag, producing a predominantly non-strange state labeled as and 
a predominantly strange state labeled as /q . Assignment of the mentioned mixed states to physical resonances is 
ambiguous because, as already discussed, there are five physical states all of which could, in principle, be candidates 
for ft and f« . 

Similarly, the isospin triplet ao can be assigned to different physical resonances - although, in this case, there are 
only two candidate states: ao(980) and ao(1450). An analogous statement holds for the scalar kaon Aq that can be 
assigned to the resonances A"q(800) or A'q(1430). A preliminary examination of the assignment of the ao field has 
been performed in Ref. [HJ where it was concluded that it most likely corresponds to the ao(1450) resonance [or, 
equivalently, ao(1450) rather than ao(980) was favored to represent a qq state]. The discussion in Ref. [l5[ was limited 
to non-strange mesons. In this work, besides the assignment of ao, we also have to consider possible assignments for 
the strange scalar field AJ; there are also two possibilities: either Aq(800) or A"q(1430). 

In the Lagrangian ([TJ there are two terms which break the original U(3)l x U(3)r [= U(3)v x ^(3)^] symmetry, 
namely the determinant term, which breaks the U(1)a symmetry, and the explicit symmetry breaking terms of 
Eqs. © and (fflh which break U{3) A , if H , A ^ and U(3) v -> SU{2) V x U{l) v , if Hg, A 8 ^ [for more details 
see, e.g. Ref. |19|]. Besides this explicit symmetry breaking the chiral symmetry is also broken spontaneously. If 
isospin symmetry is exact, only the <jm and ag scalar fields, carrying the same quantum numbers as the vacuum, can 
have nonzero vacuum expectation values (vev's) 3 . Moreover, the parameter c\ describes the axial anomaly, i.e., the 
U{1)a symmetry is explicitly broken by this term. 

As a last remark about the Lagrangian we stress that, when including a dilaton field, only two terms are modified 
[in addition, one has to add the dilatation potential, for details see Ref. [28|]: the two parameters m\ and m\ should 
be replaced by 

ml -> aG 2 , ml -> bG 2 , (8) 

where G is the dilaton/glueball field and a and b are dimensionless constants. Obviously, upon condensation of the 
glueball field, one obtains = aG\ and m\ = bG\. 

B. Tree-level masses and mixing terms 

After spontaneous symmetry breaking, the fields with nonzero vev's are shifted by their expectation values, namely, 
a N/s ~^ a N/s + 'Pn/s^ where we have introduced 4>n/s = ( a N/s)- After substituting the shifted fields into the 
Lagrangian JTJ), one obtains the tree-level masses by selecting all terms quadratic in the fields, 

C^ d = ~S i [S ij n + (m 2 s ) ij ]S j 
- \Pi [Si, □ + (m 2 P )ij] P 3 

[(-<r □ + PP)^ - rK)«] v jv 



In case of isospin breaking, also 0-3 would have a nonzero vev. 



(i 



- 2 V v {mfijk^kd* 1 ) Sj - ^Si {igifijk&kd") V jv 

+ r (gidijk<t>kdi*) Pj - -Pi {gid l]k <t> k d u ) A jv , (9) 

where {fn s )ij, (m|,),j, ('my)ij, and (rrij^ij are the scalar, pseudoscalar, vector, and axial-vector (squared) mass ma- 
trices, respectively, see their explicit expressions in Appendix El Moreover, fijj. and dijk are the antisymmetric and 
symmetric structure constants of U(3). The (squared) mass matrices are in general non-diagonal due to the mixing 
among particles sitting in the center of a given nonet, and they can be diagonalized by appropriate orthogonal trans- 
formations (for details see Appendix |A|. Besides the mixing inside the nonets there are other mixing terms, namely 
the last four terms of Eq. ©, which mix different nonets. 

In order to eliminate the latter, one performs the following shifts of the (axial-)vector fields: 

fut/s ~^ fwfS + Z m/s w hN/s d ^N/S , a>i ±,() — ► ai ±fi + Z^wa^ir^ , 
K ^±.o,o _^ R ,±.o,o + ZkWKi q, k ±o,o t K *»±fi,o _^ K *„±,o,o + z K +w K .d»K* ±firo . (10) 

These shifts produce additional kinetic terms for the pseudoscalar fields. In order to retain the canonical normalization 
for the latter, one has to introduce wavefunction renormalization constants, 

7T Znir , K Z K K , T) N / S -> Z rlN / vs T] N / s , A -> Z K *1\ p . (11) 

For the sake of simplicity we have grouped together the isotriplet states with the notation n^' , a^' and the isodoublet 
states with the notation A ± ' ^, J{*^' ' , where refers to K°. The coefficients Wi and Zi are determined in order 
to eliminate the last four mixing terms in Eq. ((§]) and to obtain the canonical normalization of the tt, fyjy, rjs, K, and 
A'q fields. After some straightforward calculation one finds the explicit expressions: 

gi4> N Vzgi&s igi(<j>N - V^4>s) gi(^N + V^4>s) nn \ 

W hN = W ai = W fls = —2 , W K * = —2 , W Kl = —2 , (12) 

rn ai '"/is *iii>K* z " t A'i 



2 m 



9l4>N 



J K 



Z J1S — 



m fl 



JlS 



2gi 



2 ^s 



Z K * 



yJim^-gl^N + V^s) 2 ' 
2rriK* 

^m^-gl^N-^s) 2 



(13) 
(14) 



It can be seen from the expressions of the wavefunction renormalization constants that they are always larger than 
one. Finally, using the explicit expressions found in Appendix [X] the tree-level (squared) masses for the different 
nonets are obtained as follows: 



l-ON 



(j>N 



J K 



ml + + 
ml+ ( Ai + 



A2 
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A2 
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V2 



4>n4>s + (Ai + A 2 ) 4> 2 s 



0i + Ai^ + ci x 



J NVS 



zi 



io + Ai^ + (A 1 +A 2 )0|+ C 



N 



- 7 2 



hps 

4>s 



n QN ,2 ,2 

h Ci <p N <p s 



vns 



are the (squared) pseudoscalar masses, while 



m o + ( Ai + -A 2 



t> 2 N 



A, 



A 2 



m o + ( Ai + — ) r N + -j=Ms + (Ai + A 2 ) 4>s 



(15) 

(16) 

(17) 

(18) 
(19) 

(20) 
(21) 
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2 2 



3 Ai 



A2 
2 



Ai</>| , 



= ml + \i(f> 2 N + 3 (Ai + A 2 ) 4>s > 



= 2A^ 



(22) 

(23) 
(24) 



are the (squared) scalar masses. The quantities m 2 NS and rn 2 TNS are mixing terms in the non-strange-strange sector. 
These mixings can be removed by orthogonal transformations, and the resulting mass eigenstates are found to be 

1 



l fa H /fo L 



))) 



n'/v 



ctn 



m ls ± \l ( m ° 



4m4. 

&NS 



K 



Moreover, the (squared) vector masses are given by 



m K * = m 1 



-(h + h 2 + h 3 )4> 2 N + 
i (gl + 2h x + h 2 ) 4>% 



1 



N 



9l) 



in 



2 



In 



us 



AT 



hi 



+ h 2 + h 3 )(j) 



25 & 



while the (squared) axial-vector meson masses are 



Ami 



h x + h 2 ) 



(25) 
(26) 

(27) 
(28) 
(29) 
(30) 



m 2 + -(2gl + hi + h 2 - h 3 ) 



26 



N 



l fis 



(gf + 2hi + h 2 



V2 



4>N4>s{h 3 - gl) + - (gl + hi + h 2 ) <j> s + 5 N + 5 



hi 



'N 



hi 



28 



s ■ 



(31) 
(32) 
(33) 
(34) 



It is interesting to note that in case of vectors and axial-vectors there are no mixings in the non-strange-strange 
sector. 



C. Parameters 



The Lagrangian (|TJ) contains 18 parameters: 

ml, ml, ci, 5 N , 5 S , gi, 52, 53, 54, 55, 56, h 0N , h ns , hi, h 2 , h 3 , Ai, A 2 . (35) 
We make the following observations regarding the model parameters: 

• The parameters hoN an< i ^os model the explicit breaking of the chiral symmetry (ESB) in the (pseudo)scalar 
sector via the term Tr[i7($ + they are uniquely determined from the mass terms of the pion, Eq. (|15[). 
and of r/s, the strange part of the r\ meson [see Eq. (|18[) ]. implying that the masses of 7? and rjs are generated 
by ESB. 

• The parameters 5n and 5s model the explicit symmetry breaking in the vector and axial-vector channels. The 
ESB arises from non- vanishing quark masses and therefore we employ the correspondence 6 n cx m 2 u d and 
5s oc m 2 s . However, in the vector-meson mass term Tr[(m^/2 + A)(i^ + RfJ] only the linear combinations 
ml/2 + S[f/s appear. Therefore, it is possible to redefine ml/2 —> ml/2 — 5n- Then, only the combination 
5s ~ <5jv appears in the mass formulas. Only this difference will be determined by the fit of the (axial-)vector 
masses. Alternately, we may set 5n — from the beginning, and determine 5s from the fit. 

• The parameters 53, 54, 55, and g§ do not influence any of the decays to be discussed in this work and are 
therefore not considered in the fit. 
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Consequently, we are left with the following 13 parameters: 

ml, m\, ci, 8 S , gi-, g 2 , h 0N , h os , hi, h 2l h 3 , Ai, A 2 . (36) 

Their large- N c dependence reads (l5| : 

ml, ml (x N° ; h N, hos oc A c 1/2 ; gi, g 2 oc N^ 1 ^ 2 ; A 2 , /i 2 , /13 oc N' 1 ; Ai,/ii oc A c 7 2 ; ci cx TV^ 3 . (37) 

In this respect the parameters 
1 . Similarly, the axial-anomaly 
model parameters (|37|) implies 
increasing N c . For this reason, 



III. FIT: RESULTS AND CONSEQUENCES 
A. Input and constraints for the fit 

Let us now turn to the fit procedure for the parameters discussed above. Our fit aims to ascertain whether (i) it 
is possible to find a fit containing masses and decay widths for (pseudo) scalar and (axial-)vector degrees of freedom 
present in our model and (ii) which physical scalar states are best described by the qq states of our model. 

All decays considered are two-particle decays, thus they can be calculated from the appropriate tree-level terms of 
the Lagrangian |T]) after applying the necessary field shifts (flO|) and orthogonal transformations (|All|) . If the decaying 
particle is denoted by A and the decay products by B and C, respectively, the decay width reads 

T A ->bc =J-^L \Ma->bc I 2 , (38) 
8irm A 

where k is the three- momentum of one of the resulting particles in the rest frame of A and Ma^bc is the transition 
matrix element (decay amplitude). Moreover, I refers to the so-called isospin factor, which is the number of sub- 
channels in a given decay channel (for instance, if B = C = K the A — > KK decay can have two sub-channels, namely 
K + K~ and K°K°, which results in 1 ~ 2). Equation (|38f will be used to calculate all decay widths entering the fit; 
details of the calculations can be found in Appendix [B] Moreover, when identical particles emerge in the final state, 
the usual symmetry factors are included. 

A brief remark is necessary regarding the errors used in the fit. The fit will contain input information from 
experimental data regarding both (axial-)vector and (pseudo)scalar states. The data are very precise for some 
of the resonances described by our model. For example, the mass of the </>(1020) resonance (our ujs state) is 
known with 0.002% accuracy. Our model does not aim to describe hadron vacuum phenomenology with this 
extreme precision. The reason is simple: already isospin-breaking effects in the physical hadron mass spectrum 
are of the order of 5% [for instance the difference between the charged and neutral pion masses, or the masses of 
the ai(1260) and the /i(1285)], but are completely neglected in our model. We therefore artificially increase the 
experimental errors to 5%, if the actual error is smaller, or we use the experimental values, if the error is larger than 5%. 

Let us now discuss the input information for our fit. We do this separately for mesons of different spin [central 
values from PDG [l[ unless otherwise stated, errors according to the above discussion]: 

• Weak decay constants: We use 

U = (92.2 ± 4.6) MeV, f K = (155.6/v^± 5.5) MeV . (39) 

The following formulas relate the decay constants to the vacuum condensates: f n = ip^/Z^ and Jk = 
(V2</>s + <f> N ) /(2Z K ). 

• Pseudoscalar mesons: We use 

m n = (138 ± 6.9) MeV, m K = (495.6 ± 24.8) MeV, 

m v = (547.9 ± 27.4) MeV, m v = (957.8 ± 47.9) MeV . (40) 

with the following notes: (i) m^ and mK represent isospin-averaged values; (ii) the relatively large error values 
come from the criterion max(5%, experimental error), discussed above. 



We recall that a non-suppressed n-meson interaction vertex scales as N c [30(. 
h\ and Ai are suppressed in the large- N c limit: in fact, they scale as -/V" 2 and not as N~ 
parameter c\ scales as iV~ 3 and not as N~ 2 . Note that the large- N c behaviour of the 
that the decay widths (the formulas for which are given in Appendix [Bj decrease with 
as already mentioned, the states in our model represent qq states. 



9 



• Vector mesons: We use 

m p = (775.5 ± 38.8) MeV, m K * = (893.8 ± 44.7) MeV, m = (1019.5 ± 51) MeV, 
T p ^„ = (149.1 ± 7.4) MeV, T K ^ Kn = (46.2 ± 2.3) MeV, T^kk = (3.54 ± 0.178) MeV , (41) 

with the following notes: (i) we use the isospin-averaged value for mx*', (ii) in case of the (j) decay we use the 
physical mass values in the kinematic factor = (m^ — 4m 2 K ) 3 / 2 due to the proximity of m^ W 2o) to the KK 
threshold (in order to eliminate phase-space effects). 

• Axial-vector mesons: We use 

m ai = (1230 ± 61.5) MeV, m /l( i 42 o) = (1426.4 ± 71.3) MeV , 
T ai ^p* = (425 ± 175) MeV, T ai ^ = (0.640 ± 0.250) MeV, r /l(142 o)^*if = (43.9 ± 2.2) MeV , (42) 

with the following notes: (i) the value of r ra] _>,,7r is not precisely known; there are experimental data suggesting 
this decay channel to be dominant for a\ |3l| and thus we estimate the possible range for r^-^ from the 
interval for the full a x decay width [= (250 - 600) MeV]; (it) according to PDG the channel /i(1420) K*K is 
dominant within the channel /i(1420) — > KKir, with the latter the overall dominant decay mode for /i(1420); 
we have assumed r / 1 (i42o)-> kk-k to be equal to the full decay width ry 1 ( 142 o) = (54.9 ±2.6) MeV and determined 
r h(ipxr>^K*K using an averaged branching ratio T fl{li20) ^ K , K /T fl{li20) ^ KI<7r = (0.8 ± 0.09) from Refs. [H] 
and 1331 . 



Isotriplet and isodoublet scalar mesons: The observables from Eqs. (|40[) - f|42[) will be used with any of the 
a -K* combinations [a (980) /A* (800), a (980)/A*(1430), a (1450)/A*(800), a (1450)/A*(1430)], where the 
data to be used are as follows: 

™a o( 980) - (980 ± 49) MeV, m oo(U50 ) = (1474 ± 74) MeV , 
m**(800) = (676 ± 40) MeV, m K * {im = (1425 ± 71) MeV , 

ra„(i450) = (265 ± 13.3) MeV, r^ (80 ow ft = (548 ± 27.4) MeV, rV» (1430) ^ = (270 ± 80) MeV , 



\M aQ[m ^ KK \ = (3590 ± 440) MeV |3j, |M ao(980W | - (3300 ± 166.5) MeV (34j . (43) 



We note the following: (i) the interpretation of Aq(800) as a particle is controversial; we will nonetheless 
include it into our hts in accordance with the conclusions of Ref. Q; (ii) as in the case of 0(1020), we will 
use the decay amplitudes rather than the decay widths for the processes ao(980) — > KK and ao(980) — > i]ir 
due to the proximity of ao(980) to the KK threshold. The reader may find it somewhat surprising that we 
are considering pairs of states above [a (1450)/A,J(1430)] 1 GeV, below [a (980)/A(J(800)] 1 GeV, as well as 
"mixed" pairs cio(980)/A'q (1430), ao(1450)/A~Q (800). The reason to consider also mixed pairs is that we want 
to avoid any kind of prejudice in the assignment of our qq states, and thus explore all possibilities. Nonetheless, 
these mixed pairs as well as the pair below l GeV will be disfavored by our analysis (see below). Mass formulas 
used in the fit are presented in Eqs. (|T5"j) - (|T5|) . (|20[) - (|23|) . and ([23]) - (|34[) . whereas formulas for the decay 
widths are given in Appendix |B| 

Moreover, the fit shall be constrained by the following conditions: 

• Z^, Z k , Z ns , Z K * > 1, due to Eqs. ^ - (H2). 

• m, lN < m ns and m aN < m as , i.e., pure non-strange states should be lighter than pure strange states. 

• uiq < 0. This is a necessary condition for the spontaneous breaking of chiral symmetry. 

• A 2 > and Ai > — A2/2. This is necessary for the potential in the Lagrangian ([1]) to be bounded from below. 

• Tn\ > 0. The reason is that otherwise (i) an instability of the vacuum in the physical p direction would occur 
[see the Lagrangian (JXJ)] and (ii) m p and m ai would become imaginary in the chiral transition, i.e., once the 
condensates vanish [see Eqs. (j2"T)) and (f3"Tj) ]. 

• mi < m p as otherwise (hi + h-2 + h^)(j) 2 v /2 + /ii0|/2 in the p mass term (|2T|) would be negative; this would 
imply that spontaneous chiral symmetry breaking decreases the p mass. This is clearly unnatural because the 
breaking of chiral symmetry generates a sizable constituent mass for the light quarks, which is expected to 
positively contribute to the meson masses. This positive contribution is a feature of all known models (such as 
the Nambu-Jona-Lasinio model and constituent quark approaches) . Indeed, in an important class of hadronic 
models [see Ref. [35j and refs. therein] the only and obviously positive contribution to the p mass is proportional 
to <t> 2 (i.e., ui\ = 0). 



10 



Before discussing the results of the fit, it is important to stress which particles have not been included in the list 
above and why. Namely, we have omitted experimental information about the scalar-isoscalar and the axial-kaon 
states. 

(i) Scalar-isoscalar states are not included in the fit because our theoretical approach is incomplete. There are 
three scalar states above 1 GeV, but at most two can be described within our model. Various approaches have shown 
(2(| that a minimal model in this energy region requires three fields, where the additional state is the scalar glueball. 
Below 1 GeV the resonance /o(500) is still poorly known and /o(980) is distorted by the nearby KK threshold. 

(ii) The axial-kaon state K\ is not included in the fit for a similar reason: the kaonic states from our axial-vector 
nonet 1 ++ mix with the kaonic states of the nonet l" 1 which is not part of our Lagrangian. (This is possible, 
because charge conjugation is not a well-defined quantum number for kaons.) Also in this case we cannot assign 
our theoretical axial kaon field to a specific resonance [in PDG there are two axial kaons with masses Ax (1270) and 
Ai(1400)[, because the mixing is expected to be large, see e.g. Ref. j36| . 

B. Results of the fit 

The experimental quantities discussed in the previous subsection do not depend on all 13 parameters of Eq. ()36l) . 
but on ttiq, Ax, mi, hi only through the two combinations 

Ci = ml + Ax (<& + 0|) , C 2 = m\ + ^ + 4) . (44) 

Moreover, instead of the parameters hoN and hos we use the condensates 4>n and <ps'- this is equivalent, as hoN, hos 
are completely determined by the masses of pion and 77s, cf. Eqs. (|15p . (|18p . Summarizing, we have the following 
eleven parameters entering the fit: 

Ci, C 2 , ci, S s , 51, g 2 , 4>n, 4>s, h 2l h 3 , A 2 . (45) 

We perform our fit using the experimental values for the 17 quantities given in Eqs. (|39j) - (|4"2"|) . In addition, we use 
the experimental values for the pairs a (980) / A*(800), a (980) / #£(1430), a (1450) / A*(800), a (1450) / A*(1430) 
given in Eq. (|43|) . These are four additional experimental quantities [except in the case of ao(980), where we use the 
two values for the decay amplitudes instead of one value for the decay width]. In total, we therefore fit 21 (22) 
experimental quantities to the eleven parameters given in Eq. (|45j) . The corresponding values of x 2 are listed in 
Tablefl] We see that the combination ao(1450)/A'o(1430) gives the best value for \ 2 ■ In fact, with a x 2 e d of 1.23 this 
fit is remarkably good, meaning that all physical quantities are reproduced within an average error of 5%. In the 
framework of an effective model for the strong interaction, which neglects isospin-breaking effects and describes all 
mesonic resonances up to about 1.8 GeV, we perceive this to be a remarkable achievement. 

Note that also the fit with the pair ao(980) / Aq(1430) has a Xred °f the same order of magnitude as the one 
with the pair ao(1450) / Aq(1430). Nevertheless, we will disregard this fit for two reasons. Firstly, the Xred f° r the 
pair ao(980) / Aq(1430) is, although rather small, still larger by a factor of about two than the Xred f° r the pair 
ao(1450) / ATq(1430). Secondly, and more importantly, the fit with the pair ao(980) / Kq (1430) produces a scalar-kaon 
mass of mx* = 1146 MeV, which cannot be assigned to any physical resonance as it is much larger than m K *( S00 ) 
and much smaller than tuk* (1430I: Note that this problem is also present in NJL models with mixing between scalar 
mesons below and above 1 GeV |37| . 



Pair 


x 2 


Xred 


ao (1450)/ Kq (1430) 


12.33 


1.23 


o (980)/ifo(800) 


129.36 


11.76 


a (980) /Kq (1430) 


22.00 


2.00 


a ( 1450) /.Kg (800) 


242.27 


24.23 



TABLE I: Isotriplet and isodoublet scalar pairs and the corresponding values of the total \ 2 an d the reduced Xred = X 2 /Ndof, 
where Ndof is the difference between the number of experimental quantities and the number of fit parameters (10 for the first 
and fourth rows and 11 for the second and third rows). 

The detailed comparison of theory with data for the pair ao(1450) / Aq(1430) is presented in Table ITT1 where the 
theoretical errors are also shown. Errors for the model parameters (Spi) are calculated as the inverse square roots of 
the eigenvalues of the Hessian matrix obtained from x 2 (Pj ) > where pj denotes elements of the parameter set (|4"5)) and 
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i = 1, 2, . . . , 11. Theoretical errors AOi for each observable Oi (mass, decay width) can be calculated according to 
the following formula: 



AO,- = 



\ 




S P:j 



(46) 



at fit value of Oi 



The remarkable agreement of our results with experimental data in the (pseudo)scalar and (axial-)vector sectors can 
now be explicitly seen for the various quantities shown in Table [n] We thus conclude that the quark-antiquark state 
a,Q and Kq should be assigned to the resonances ao(1450) and ^(1450). Moreover, also the axial-vector resonances 
are well described by the fit: we therefore also interpret them as predominantly quark-antiquark states. In Table ITTT1 
we finally present the values of the parameters of the model together with their errors as obtained from our best fit. 



Observable 


Fit [MeV] 


Experiment [MeV] 


u 


96.3 ±0.7 


92.2 ±4.6 


fK 


106.9 ±0.6 


110.4 ±5.5 


m-n 


141.0 ± 5.8 


137.3 ±6.9 


m K 


485.6 ± 3.0 


495.6 ± 24.8 


m v 


509.4 ±3.0 


547.9 ± 27.4 


«v 


962.5 ±5.6 


957.8 ±47.9 


m p 


783.1 ± 7.0 


775.5 ± 38.8 


m K * 


885.1 ±6.3 


893.8 ±44.7 




975.1 ±6.4 


1019.5 ±51.0 


m ai 


1186 ±6 


1230 ± 62 


m /l(1420) 


1372.5 ± 5.3 


1426.4 ± 71.3 




1363 ± 1 


1474 ± 74 


m *s 


1450 ± 1 


1425 ± 71 


Tp— >7T7r 


160.9 ±4.4 


149.1 ± 7.4 




44.6 ± 1.9 


46.2 ±2.3 


f*-tKK 


3.34 ±0.14 


3.54 ±0.18 


I^ai — y p-7r 


549 ± 43 


425 ± 175 


Tai — >7T"y 


0.66 ±0.01 


0.64 ± 0.25 


r/l(1420)->K*if 


44.6 ± 39.9 


43.9 ±2.2 




266 ± 12 


265 ± 13 




285 ± 12 


270 ± 80 



TABLE II: Best-fit results for masses and decay widths compared with experiment. 



Parameter 


Value 


Ci [GeV 2 ] 


-0.9183 ± 0.0006 


C 2 [GeV 2 ] 


0.4135 ±0.0147 


ci [GeV -2 ] 


450.5420 ± 7.0339 


5s [GeV 2 ] 


0.1511 ±0.0038 


gi 


5.8433 ± 0.0176 


92 


3.0250 ± 0.2329 


<Pn [GeV] 


0.1646 ± 0.0001 


0s [GeV] 


0.1262 ±0.0001 


hi 


9.8796 ± 0.6627 


h 3 


4.8667 ± 0.0864 


A 2 


68.2972 ± 0.0435 



TABLE III: Parameters and their errors. 
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C. Consequences of the fit 

1. Scalar-isoscalar mesons 

We now turn to the scalar-isoscalar mesons. We shall discuss four different aspects: 

(a) Results in the large-N c limit (X\ = h\ = 0): 

From the fit in the previous section we cannot immediately obtain the masses of the scalar-isoscalar states of the 
model because their masses do not depend solely on the combination C\ = TOq + Ai {<$> 2 N + <^|), but separately 
on the parameters ?tt,q and Ax- Similarly, the decay rates of the scalar-isoscalar mesons do not depend solely on 
the combination C% = m\ + 4j- (<p 2 N + but separately on the parameters vti\ and hi. 

Interestingly, the parameters Ai and hi are large- N c suppressed. Setting them to zero we obtain a prediction for 
the masses and decay widths of the scalar-isoscalar states in the large- limit. Moreover, when Ai = hi = 0, 
there is also no mixing between fin and ss states: f$ = ajv = fin = (uu + dd)/y/2 and /q = as = ss. Their 
masses read: 

m f L = 1362.7 MeV, m f H = 1531.7 MeV . (47) 

These masses clearly lie above 1 GeV. Comparing these values to the experimental values for the three isoscalar 
states above 1 GeV [l| 

m/ (i370) = (1350 ± 150) MeV, m /o{1500) = (1505 ± 75) MeV, m /o(mo) = (1720 ±86) MeV , (48) 

we see that the mass of ffr is well compatible with that of the resonance /o(1370). [Note that m/ (i37o) in Eq. 
P8"]) emulates the PDG mass interval (1200-1500) MeV.] The mass of / ff appears to be within that of / (1500). 
Remember, though, that we artificially enlarged the experimental error to 5%; the actual error is only 6 MeV, 
and then the mass of /q would be (just) outside the experimental error. Considering the decays of /q we shall 
provide evidence that an assignment of /g to /o(1710) is also possible (and even rather likely). 

For Ai = hi = the decay rates of f$' H into tttt and KK read 

r f L ^ = 520 MeV, r f L _^ KK = 129 MeV, (49) 
T = MeV, T f H^ KK = 422 MeV. (50) 

For the experimental values we quote the PDG values [l|: 

r/o(1370)->7T7r = (250 ± 100) MeV, ^f (l370)^KK < r/ (1370)->-7r7r j 

r/ o( iBoo)-Hnr = (38 ± 5) MeV, r /oC i B oo)-,jcjc = (9.4 ±2.3) MeV , 

r /o( moW = (29.3 ±6.5) MeV, r /o(m0 )^ K = (71.4 ±29.1) MeV . (51) 

Note that for the resonance /o(1370) no branching ratios into irn and KK are reported in Ref. [l|. The value for 
r/o(i370)->7T7r is our estimate from PDG and the review Q, whereas ^f n37o)-^KK is our estimate from results 
presented in Refs. P. [38|. 

Our results for f$ are in agreement with the experimental decay widths of /o(1370). Our theoretical value for 
f$ — > KK turns out to be too large, while /q — > irir vanishes, because /q = 175 is a pure ss state. Nevertheless, 
our model predicts the existence of a scalar-isoscalar state which decays predominantly into kaons; this is indeed 
the decay pattern shown by /o(1710). For these reasons we identify our state f$ as (predominantly) /o(1710). 
This identification will be further tested when Ai 7^ 0, hi 7^ below. 

It should be stressed that a quantitative study of the scalar-isoscalar system cannot be performed at present 
because our model contains only two states, while three resonances appear (this is also the reason why we did 
not include the scalar-isoscalar states into the fit). As many studies confirm [26[, the mixing with the scalar 
glueball can be sizable: a reliable analysis of the scalar-isoscalar states can only be performed when taking the 
scalar glueball into account. The preliminary results in the Nf = 2 sector [28| have indeed shown that the 
glueball and quarkonia degrees of freedom interact strongly and that the decay patterns are sizably influenced. 

(b) Results for Ai ^ , hi + 0: 

A nonvanishing value of the parameter Ai induces a mixing of the pure states ojq and erg, see Eq. (|2"4")h Our 
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fit determines only the value of the linear combination C\ rather than the value of Ai, see discussion before 
Eq. (|4"5)) . Nonetheless, a range of values for Ai can be estimated using the value of C\ from the fit and the 
condition tUq < 0. Consequently, we also obtain a range of values for our isoscalar masses rrijL.H and isoscalar 
decay widths. The masses vary in the following intervals: 415 MeV < m^z, < 1460 MeV and 1480 MeV 
< rrifH < 1981 MeV. Considering these mass values, /q may correspond to either /o(500), /o(980), or /o(1370) 

and /o* may correspond to either /o(1500) or /o(1710) [![. Therefore, a mere calculation of scalar masses does 
not allow us to assign the scalar states f$ and /q to physical resonances. In order to resolve this ambiguity, 
we will calculate various decay widths of the states and f$ and compare them to data [l[ . 




The decay width Y f L, is consistent with the experimental range (250 ± 100) MeV for the /o(1370) — > irir 

J 

decay width, if 1000 MeV < m^L < 1460 MeV. Other assignments can be excluded: (?) our f$ state cannot be 
assigned to /o(500) as T f L^„ < 20 MeV in the (new) PDG mass range (400 - 550) MeV for / (500). Therefore, 
this strongly disfavors /o(500) as a qq state, (ii) TjL_ ¥mt varies between approximately 130 MeV and 150 MeV 
in the mass interval of the /o(980) resonance, i.e., the interval between 960 MeV and 980 MeV. The PDG result 
for the full decay width of the Jo (980) resonance is between 40 MeV and 100 MeV, with the mr channel being 
dominant jl|, i.e., about a factor two smaller than our theoretical values. Consequently, the assignment of Jq 
to /o(980) is also disfavored by our analysis. Therefore, (although there could still be mixing with other states) 
we assign our state to /o(1370), thus supporting the interpretation of the latter state as a quarkonium. 
Although experimental data is not conclusive, we mention that the decay width Tjl^ kk shown in Fig. [2] is 
consistent with the data. 

In the mass interval 1500 MeV < f$ < 1560 MeV, our results for the decay width Tjh^^ are consistent with 
experimental values for both /o(1500) and /o(1710), see Eq. (|5ip . However, the KK decay widths present us 
with a different conclusion: as already indicated in the case Ai = hi = 0, our /q state decays much more 
abundantly into kaons than into pions, and the experimental data suggest only one physical resonance with 
the same feature: /o(1710). The other resonance /o(1500) decays preferably into pions. Thus (and although 
the mass value for is too small in the range where the decay width into pions agrees with the experimental 
value) we assign our f$ state to /o(1710) and inspect in the following whether this assignment is justified by 
other data. 

We observe from Fig. [2] that r fH^ KK rises rapidly and remains above the PDG result T f (nity_>KK = (71.4 ± 
29.1) MeV in the entire mass interval TUfH > 1500 MeV. Although the absolute value of the decay width is rather 
large, several ratios of decay widths can be described correctly by our fit, most notably the irir/KK decay ratio 
presented in Fig. [3] 

Let us first discuss results for r fL^. KK /T fL_ yvK (left panel of Fig.[3|). We observe that the ratio varies between 
0.49 for m f L = 1152 MeV and for m f L = 1444 MeV. Experimental data regarding this ratio for / (1370) 
are unfortunately inconclusive, ^f (i37o)^-K K /r/ (i37o)-Mr7r = 0-08 ± 0.08 is quoted by the BESII Collaboration 
(39l | and Fj ( 137 o)_ ! .iv"/<"/r/ (i37o)->7r7r = 0.91 ±0.20 is the value given by the OBELIX Collaboration [4(|; for the 
mass interval shown in Fig. [31 our result is most consistent with the value of the WA102 Collaboration [4l| . 
rj ( 13 7o)_j.xA" /r/ (i37o)->-7r7r = 0.46 ±0.15±0.11. The ambiguities in the experimental value of this ratio do not 
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FIG. 2: Y jl^ kk and Y ^h^ kk as functions of m^L and m^n, respectively. 




allow us to constrain our parameters, although our f$ state is compatible with the /o(1370) data also for this 
particular case. 

Let us now discuss results for r fH^. nn /T fH^ KK (right panel of Fig. [3]). As already mentioned, the analysis 
of the two-pion decay did not allow for a definitive assignment of our f$ state, as it could correspond ei- 
ther to / (1500) or to / (1710). The PDG data suggest the ratio Y fo(15m) _, KK /Y fo(15m) ^^ = 0.246 ± 0.026 
or ry ( 150 o)-». W 7r/r^ (i5oo)->.jfjf = 4.065 ± 0.430, respectively, whereas, for / (1710), we use the WA102 ratio 
r/ (mo)^7r7r/r/ (i7io)-> ifi? = 0.2 ± 0.06 [4l| . The latter is only marginally consistent with the one preferred by 
the PDG lr™f noHTI[ /r™ G noHJf}f = 0.4lt°;^], originally published by the BESII Collaboration [42j, that 
suffers from a large background in the Tr + Tr~ channel (approximately 50%) and is therefore omitted from our con- 
siderations. We observe from the right panel of Fig.|3]that our ratio r jh _^ nv /T jh _^xk never corresponds to the 
one experimentally determined for /q(1500). Although the ratio shows a strong increase near the left border of the 
mass interval, decreasing the mass beyond this border would violate the constraint m§ < 0. Thus, the experimen- 
tal value is out of reach. Conversely, our ratio describes exactly the value ry ( 1710 )_ i . 7I . 7r /rj ( 1710 )_^.^^ = 0.2±0.06 
if = 1502^3 MeV. [We discard the second possibility rrifB = I6HI23 MeV as the decay width Y ^h^ kk 
would exceed 600 MeV, see Fig. [21 i.e., an order of magnitude larger than the experimental value, see Eq. ([51]) .] 
This implies Ai = —4.1 qp 0.7 - the parameter Ai remains close to the large- N c limit and is much smaller than 
A 2 , see Table HTU 

The contribution of the pure-strange state as to /q is then approximately 96%, as can be calculated from Eqs. 
(|A12[) - (jAl6[) . Thus our predominantly strange state /q describes the irw/KK ratio of /o(1710), and not 
that of /o(1500), although the correct description requires that TUfs corresponds to the mass value of /o(1500). 

The correct description of the decay ratio is an indication that our assignment of /q to /o(1710) is the correct 
one, whereas the fact that TrifH is smaller than f7i/ (mo) indicates the necessity to include a third isoscalar 
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degree of freedom into our model: the glueball. As concluded from the Nf = 2 version of our model [28fl, the 
/o(1500) is predominantly a glueball state and thus the inclusion of the glueball state is expected to induce a 
level repulsion in the masses. This may shift the currently too small mass value of the predominantly strange 
quarkonium from approximately 1.5 GeV to 1.7 GeV, where it is experimentally found. Nonetheless, our results 
clearly demonstrate that scalar quarkonia are found above, rather than below, 1 GeV. 

Until now we have only considered the case where one of our large- N c suppressed parameters (Ax) is nonzero. 
We have also investigated the influence of non-vanishing values for the other large- N c suppressed parameter, 
hi, on the decay widths lY^-Mm and r jl.h ^ kk [see Eqs. (|B24|) . (jB25[) , (|B30[) , and (|B31[) ]. Consistency with 
the large- N c deliberations requires us to keep hi smaller than, or in the vicinity of, hi and /13, see Table UTTl We 
again observe that our ratio T jh^^^/T jh^ kk never corresponds to the one of the / (1500) resonance whereas 
the ratio r / (mo)-Hr7r /F f (mo)^KK = 0.2 ± 0.06 is correctly described if hi ~ —9. In this case, nrifH rises to 
approximately 1540 MeV and is thus outside the PDG result m,y ( 15 oo) = (1505 ± 6) MeV but still too small 
when compared to m./ (mo) = (1720 ±6) MeV. Thus, it is still necessary to include a glueball degree of freedom 
into our model. The assignment of f$ to /o(1710) is further justified considering decays into r\ and 7/ mesons. 

(c) r\r\ and r\vl decay channels for the scalar-isoscalar mesons: 

PDG data suggest the following values of r\r\ decay widths for /o(1500) and /o(1710) [l| 

r /o(1500Hm = (5.56 ±1.34) MeV, r /o(mo) ^ w = (38.6 ± 18.8) MeV . (52) 

Our analysis always yields T fH^ m > 20 MeV; there is consequently no value of at which T jh^^ would 
be compatible with r / (i5oo)->rjrj- However, tu^h = 1502^ MeV [obtained requiring r fH_ yw „/T fH^ KK = 
r /o(i7io)^7 r7 r/r/ (i7io)->/<'^ J i-e., Ai = -4.1 =F 0.7 and hi = 0] yields 

r /f _ MW = 49.61^3 MeV, (53) 

(where the errors arise from the uncertainty in Ai only), i.e., in agreement with the experimental value quoted 
in Eq. (|52j) . The same parameter set also yields 

r^.^- 33 MeV; (54) 

purely a prediction, as ^ f (i37o)-i. nn has not yet been measured. 
The choice Ai = —4.1 =p 0.7, hi = also yields 

T fk^-nnl^ fk^KK = °-194±a°o3 , T f H JT f H^ KK = 0.204t";™2 , (55) 



/o fo^KK — ' L fff^vvl f"^ KK ~ -0.002 i 

fo-trrn/FfL^Tm = 0.0871^006 ' ^■fo^rm/^-fa^'KT' = ^-^+0.42 ' 



whereas experimental data read 

r f (i500)^rin/ T Mi500)^KK = 0.59 ±0.12 Q , (57) 

r /o (i5oo)^ w /r /o (i5oo)^ - 0.145 ±0.027 i, (58) 

r /»(mo)^/r/„(ino)^ = 0.48 ±0.15 [H , (59) 

r /o (mo)W r /o(mo)^ = 0.461°;™ , (60) 

r/ (mo)^/r /o (mo)^ = 2.40 ±1.04 [MIH . (61) 

In all cases, the /o(1500) data are off by several standard deviations from our theoretical results, while there is 
good agreement with the results for /o(1710). 

We have also considered decays involving the n' meson. Since the threshold for 777/ decays is at approximately 
1.5 GeV, it suffices to consider decays of f$ only. The value Ai=— 4.1=f0.7 yields 

T f n^ m , = 12.7±i:i MeV , (62) 

r /rf->W'/ r /o ff ^^r = °- 26 +0.04 > ( 63 ) 

r /o ff ^W'/ r / ( f ->W = °- 26 io.05 > ( 64 ) 

TfM^rf/TfB^KK = 0-05 ± 0.01 , (65) 
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whereas experimental data read 

iMisooW =(2.1±1.0)MeV, (66) 

r/ (i5oo)^ W '/ r /o(i50o)^7T7 r = 0.055 ± 0.024 . (67) 

F/ o (1500)->?7?7'/r/ o (1500)->-)7'7 = 0-38 ± 0.16 . (68) 

The decay ratio Y ^^qq^,^ /Y f ^ 1500 ^ KK is unknown; there are also no data for the 7777' decay channel of 
/o(1710). Still we observe that neither T jh^^, nor YfH^^i /Y fH^ vn describe the corresponding experimental 
results for / (1500). Indeed the only piece of experimental data regarding / (1500) that is described by our fit 
results is the one for the rjrj' /rjr] decay ratio. Nonetheless, all the other results regarding decay ratios obtained by 
our analysis clearly demonstrate the correspondence of our predominantly strange state f$ to /o(1710); fy was 
found to correspond to /o(1370) already in the discussion of the 7T7t decay channel. Consequently, we conclude 
that /o(1370) and /o(1710) are favoured as scalar qq states. However, we also stress again that the mass of our 
remains too low when compared to TO/ (mo) due to a missing scalar- glueball state expected to shift rrifH 
to TO/ (i7io) by level repulsion. 

2. Mixing in the pseudoscalar-isoscalar sector 

Our Lagrangian {1} implements the mixing of two pure pseudoscalar isosinglet states, t]n = (uu + dd)/y/2 and 
r/s = ss. The mixing term is presented in Eq. (|19l) . The mass terms for 77 jv (|17|) and rjs (JTSj) are determined by our 
fit parameters presented in Table ITlTl The same parameters also determine the 77-77' mixing angle as [see Eqs. (|A12[) - 

(SEN 

6 V = -44.6° . (69) 

The result is close to ideal mixing. 

3. The axial-vector kaon state K\ 

The K\ state has not been assigned to a physical resonance and included into our fit, because the PDG listing 
suggest two distinct assignment candidates: i4Ti(1270) and ifi(1400), expected to mix (36|. Therefore, mass and decay 
widths are left as predictions of this work. 

There are three decay widths of the K\ state that can be calculated within our model at tree level: K\ — > K*ir, 
pK, and u>K. They account for approximately 70% of the full ifi(1270) decay width and almost 100% of the full 
A'i(1400) decay width [l|. Using the parameter values stated in Table [TTTI it is possible to calculate the mass, Eq. 
(|32| . as well as the decay width [via the generic decay- width formula (|B46|) [ of our K\ state. We obtain the following 
results: 

m Kl = 1282 MeV , T Kl ^K*n = 205MeV , Y Kl ^ pK = 44 MeV , T Kl ^ uK = 15 MeV . (70) 

The mass is within 2a of 771^(1270) = (1272 ± 7) MeV and thus rather close to the experimental result. However, 
the value of the full decay width is 264 MeV, while the PDG data read IV^mo) = (90 ± 20) MeV and Y^Iiaoo) = 
(174 ± 13) MeV. Our result is therefore approximately three times too large when compared to the data for if 1 (1270) 
and approximately 50% too large when compared to the data for i^i(1400); errors have been omitted from the 
calculation. These results demonstrate the necessity to include a pseudovector I(J PC ) = 1(1^ ) nonet into our 
model and implement its mixing with the already present axial- vector nonet. A possible mixing term between an 
axial-vector nonet A\ and a pseudovector nonet B\ reads 

Tr(A[A1,B lfl ]) , (71) 

with A from Eq. ([7]). Various studies have indeed found the mixing to be non- negligible (3(|; see also Ref. (45| . 

4- Branching ratios o/ao(1450) 

As a consequence of our fit we can determine the branching ratio of the resonance ao(1450) into KK, nr), and 7777'. 
We obtain the following values [see Eqs. ([BT]) . (|B2]> , and (|B11[) ]: 

Y a „ = (115.4 ± 6.2) MeV, r Qo ^ = (21.5 ± 1.4) MeV, Y ao ^ KK - (128.8 ± 3.9) MeV. (72) 
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This leads to the following branching ratios: 



0.19 ±0.02, T °°^ KK =1.12 ±0.07, (73) 



r r 



which should be compared with the experimental results [1| 



r Qo( 1450)^^ = q 35 ± q lg i r ..(1450)^ = Q gg ± Q 23 _ (74) 
r ao (1450)^i77r r ao ( 1 45o)^ ?)7r 



Our results are, within errors, consistent with the data. 



5. Contributions to the mass of the p meson 



The mass of the p meson consists of three terms, 

m 2 p = m\ + \{hi+ h 2 + h 3 )tf N + -±<f> 2 s . (75) 

The first term m\ is generated by the condensation of the dilaton field G, m\ cx G\. The second term is proportional 
to the chiral condensate 4> 2 N . In the large- ./V c limit the parameter h\ = and we can determine the terms from the 
results of the fit: mi = 0.643 GeV and a/ (h 2 + h 3 )/2<fijy = 0.447 GeV. It turns out that the glueball-driven term is 
dominant. 

It is interesting to note that models based on QCD sum rules (4|| or Brown- Rho scaling [4?} predict that m p oc <Pn, 
where k = 2 in the first and k = 3 in the latter case. The order parameter for chiral symmetry breaking c/)jv decreases 
as a function of nuclear matter density n: for small densities it is well known (see, e.g. Ref. [46j) that 

_M^«l-0.3^, 

4>N I vacuum no 

so these models predict a substantial decrease of m p already at nuclear matter saturation density, no- In this respect, 
it is interesting to evaluate the value of the p mass at no according to Eq. (|75|) . We assume that the glueball-driven 
term m\ does not vary. This can be motivated by considering that the glueball is massive and a substantial decrease 
in its mass may occur only at much higher density. Assuming that the linear dependence on the density of the 
condensate holds u p to saturation density, we obtain that m p decreases by about 70 MeV compared to 90 — 130 MeV 
predicted by Refs. (MSI- 

Obviously, these considerations are only qualitative but could represent the starting point of interesting studies of 
vector mesons in the medium, which is an important aim of several experimental investigations [48[. To this end, 
one has to calculate the behavior of the dilaton and chiral condensates in the medium within the same theoretical 
framework. 



IV. CONCLUSIONS 



We have presented a linear sigma model with three quark flavors. The model implements the symmetries of QCD, 
the discrete CPT symmetry, the global chiral U(Nf)L x U(Nf)n symmetry, and the breaking mechanisms of the 
latter: spontaneous (due to the chiral condensate), explicit (due to non- vanishing quark masses) as well as at the 
quantum level [the U(1)a anomaly]. 

The model has been used to describe meson states up to energies of ~ 1.7 GeV. This energy region exhibits numerous 
resonances, related by scattering reactions and decays. For this reason, a realistic model of QCD degrees of freedom 
in the mentioned energy region should describe as many of the resonances as possible. Thus, we have constructed a 
linear sigma model that contains scalar (two isoscalars, /q and /q , as well as an isotriplet, ao, and two isodoublets, 
Kq), pseudoscalar (tt, K, ij, ij'), vector (p, w, K*, (j>), and axial- vector [oi, K\, /i(1285), /i(1420)] degrees of freedom. 
We have thus constructed a single model that contains (pseudo)scalars and (axial-)vectors both in the non-strange 
and strange channels. To our knowledge, this is the first time that such a comprehensive approach has been presented. 

The model, dubbed extended Linear Sigma Model (eLSM), has allowed us to study the overall phenomenology of 
mesons and, in particular, to explore the nature of scalar and axial-vector resonances. In order to test our model we 
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have performed a global fit to 21 experimental quantities involving both the (pseudo) scalar and the (axial-)vector 
masses and decays. Due to mixing with the scalar glueball (not included here), we did not include the scalar- isoscalar 
resonances in the fit. Similarly, we have omitted the axial- vector resonance K±, due to the fact that in reality a large 
mixing of two kaonic fields from the 1 ++ and l" 1 nonets takes place. 

One of the central questions of our discussion has been the assignment of the scalar states: to this end we have 
tested the possible scenarios for the isotriplet and isodoublet scalar states by assigning our scalar fields ao to ao(980) 
or ao(1450) and Kq to Kq(800) or A'q(1430). The outcome is univocal: the global fit works well only if the states 
ao(1450) and Aq(1430) are interpreted as quark-antiquark states. On the contrary, the other combinations deliver 
large and unacceptable values of x 2 , see Table [I] We thus conclude that the scalar 1 = 1 and 1 = 1/2 states lie above 
1 GeV and have to be identified with the resonances ao(1450) and Aq(1430). Moreover, the overall phenomenology 
described by the fit is very good, see Table [TT1 It is then possible to describe properly many different mesonic masses 
and decays within a unified treatment based on the symmetries of QCD. It should be stressed that the inclusion of 
the (axial-)vector mesons has a crucial impact on our results and represents the most important new ingredient of 
our approach. The good agreement with data also shows that the axial-vector mesons can be interpreted, just as 
their vector chiral partners, as quark-antiquark states. 

We have then studied the consequences of our fit. We have primarily concentrated on the scalar- isoscalar sector 
which was not included in the fit. In the limit N c — > oo it is possible to make clear predictions for the two states 
Jq and /q . Their masses lie above 1 GeV and their decay patterns have led us to identify /q with (predominantly) 
/o(1370) and f$ with (predominantly) /o(1710). The theoretical decay rates of f$ are in agreement with experiment; 

decays only into kaons, but turns out to be too wide. At a qualitative level, this large-iV c result clearly shows 
that also the scalar-isoscalar quark-antiquark states lie above 1 GeV. The overall situation in the scalar-isoscalar 
sector can be slightly improved when studying the influence of large- N c suppressed terms. However, in this system 
the inclusion of an additional glueball state would be necessary for a full quantitative study. Mixing phenomena 
are known to be large here and aff ect b oth the masses and decays. The study of a genuine three-state system is 
mandatory: starting from er/v = \/T/2(uu + dd), G = gg, and as = ss to describe properly /o(1370), /o(1500), 
and /o(1710). Additionally, we have studied other consequences of the fit, such as predictions for the ijij channel 
of scalar-isoscalar states, the rj-rj' mixing, the prediction for the axial-vector kaonic state, and for the ao branching 
ratios. Finally, we have also discussed the origin of the mass of the p meson and, at a qualitative level, its possible 
value at nuclear matter density. 

There is, however, one question that remains open. Firstly, interpreting resonances above 1 GeV as qq states 
leads to questions about the nature of /o(500), ao(980), /o(980), and A"q(800). Their presence is necessary for the 
correct description of mr scattering lengths [see, e.g. Ref. [15f]. There are two pos sibilities: (z) they can arise as 
(quasi-)molecular states of mr or KK, respectively, or (u) as tetraquark states [49]. The question is whether the 
attraction is large enough so that these states are bound, or whether the attraction is not sufficient so that they 
are (unstable) resonances in the continuum. This question and deciding whether possibility (i) or (ii) is realized in 
nature represent interesting starting points for future studies, along the lines of the Bethe-Salpeter ap pro ach of Ref. 
[50l | . the lattice approach of Ref. (Hl|, or even holographic approaches such as those presented in Ref. [52]. In either 
case, one may include the low-lying isoscalar states as interpolating fields in our Lagrangian, such as has been done 
in Ref. |H|. 

The present model can be also studied in the baryonic sector, see Ref. pjj] for the two-flavor case. The very same 
ideas of chiral symmetry and dilatation invariance can be applied to the baryonic sector as well. The extension of 
the model to three flavors in the baryonic sector would surely represent an interesting problem: a multitude of data 
on decays and masses is available to make a precise test of our approach. 

Additionally, restoration of chiral symmetry at nonzero temperature and density is one of the fundamental questions 
of modern hadron physics, see, e.g. Refs. pj| [HI, in which the two- flavor version of this model has been studied at 
nonzero density, or Ref. |55| for a holographic approach to the exploration of nonzero chemical potential. Linear 
sigma models are appropriate to study chiral symmetry restoration because they contain from the onset not only 
(pseudo)scalar and (axial-)vector mesons but also their chiral partners; mass degeneration of chiral partners represents 
a signal for the chiral transition. Therefore, we also plan to apply the model to study chiral symmetry restoration at 
nonzero temperatures and densities. 
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Appendix A: Tree-level masses 



After spontaneous symmetry breaking (see Sec. Ill Bj) . from the quadratic terms of the Lagrangian the (squared) 
mass matrices for the different fields are read off as 



( m s)ij = S ip {mlS pq + iF pq i m cj)i(t) m )S q j, 

(rrip)ij = Si P (mlS pq — 9ciG p i m 4>i<f) m G q i> m i<f)i>(j) m i + AH pq j m <t>i<l> m )S q j, 

{lTly)ij = Sipim^Spq + J r pq lm4'l ( t > m) S ' q j , 
(m 2 A )ij = Sip(m% g + J pq lm<t>l<l>m)S q j, 

where S q j is the (9 x 9) transformation matrix to the non strange - strange base, 



/a/2 0... 1 \ 



5' 



V3 



1 



\ 1 ... -V2 J 



and the F, H, G, J, J' coefficient tensors are 

Fijki 



-J (SijSkl + SikSjl + SuSjk) + (dijmdklm + diUmdjlrn + dil m djkm) ■ 



Gijk 



SjjSkl -\- (dij m dklm ~t~ fikmfjlm ~t~ filmfjkrn) ■ 
4 o 



9 3 
dijk + -^doooSioSjoSko — ^ wodjko + Sjodiko + Skodijo) 



J, 



ijkl 



V 

,J ijkl 



9l fikmfjlm ^ $ij$kl ~\~ ^ ^ijm^klm ~t~~ ^ {^ikm^jlm H~ &ilm&jkm fikmfjlm filmfjkm) ; 

9 > Mm + £ M « + h i^d klm - h i ( dikmdjlm + damdjkm - Mm - f M 



(Al) 
(A2) 
(A3) 
(A4) 



(A5) 



(A6) 
(A7) 
(A8) 
(A9) 
(A10) 



As can be seen from Eqs. (fl"T)) - (fT9"|) . and Eqs. (|2"2l (|2"4")h there is a mixing in the pseudoscalar and scalar N — S sectors, 
which can be resolved by the following two-dimensional orthogonal transformations 



O 



cos ( 



rj/<r 



rj/tj 



sm( 



rj/o- 



sin^yv cosO v / a 



(All) 



where 9 v / a are the pseudoscalar and scalar mixing angles. The (squared) mass eigenvalues can be written with the 
help of the mixing angles as 



to 2 i„ = mt i„ cos 2 Q„i„ + nii ,„ s'm20 r ,/ rr + m„ /„ sin 2 0„/„, 

VN/&N r ll' J VNS/&NS Vf 7 »)/< J ' 



(A12) 



l m/<T 2 ~ m VN/a N Sln V-n/a - m r )NS /<y N 



sin 26 



n/a + m nsla s 



cos 2 e. 



■q/ai 



where, if we require that w 2 ^^ > m 2 ^^, the mixing angles are given by 



-sign(m 2 ;„ 



1 

71 



m? , — m? , 



(m 



, — m 
»)n/o"jv vs 



, ) 2 + 4m 4 , 

las' riNS/a 



(A13) 
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cos6>„ /CT = ^ 



? i 
to i — Tci i 

VS/as riN/afr 



\ < /(to 2 , — to 2 , ) 2 +4m 4 , 

\ V r IN/o- N Vs/a-s' riNs/a 



(A14) 



It can be seen that, if m 2 , 
these expression into Eq. (|A12[) it is found that 



> 0, then — 7r/2 < 6 v / a < and, if \ 



VNS/&N! 



< 0, then < 6 v / a < tt/2. Substituting 



v 2 =H 

J = 1 r 



W 

Vn/<?n 



„ + m 2 ,„ — »/(m 2 , — m 2 , ) 2 +4m 4 , 



I ,„ +J(m 2 , -m 2 . ) 2 +4m 4 

ris/<?S V V i?N/ffJV Vs/a- s J V 



Vns/ "NS 



(A15) 
(A16) 



identified as to,,/, and similarly for a. 



,2 



is fulfilled. In this way we know that m V2 must be 



Appendix B: Decay widths 



In this section we show the explicit form of some of the most relevant decay widths calculated from our model at 
tree-level using Eqs. ([38]) . The formulas below are organized according to the type of the decaying particle. 



1. Scalar-meson decay widths 

At first we are considering the decays of the scalar isotriplet ao, scalar kaon Kq, and scalar isosinglets / . The ao 
state has three relevant decay channels, namely into tot, rj'ir, and KK, the first two of which are strongly connected 
due to the mixing between 77 and 7/. The ao — > r)ir and ao — > rj'ir decay widths 4 , considering that in both cases 1=1, 
read 



f 



8to Qo 7t 



8to Qo 7t 



(™« - m l ~ m l? ~ 4to 2 TO 2 



1 1/2 



4to 4 



("<, " "V - m i) - 4m 2 ,mi 



\M 



4< 



1/2 



|A4 ao — >ri'n\ 1 



(Bl) 
(B2) 



with the following transition matrix elements, 



M ao -> m = cos6 w M ao ^. m A m v) + Sm9 7T Mao->r,s^( r n r i) , 
Ma -> v 'iT = COS w M ao ^.r, 3 A m v') ~ sin ^•W«o->W» i m v' ) 



(B3) 
(B4) 



where 



•fi^aQ— lT}NTr\JYv) -A ao1 -j N7r B aQ7 j N7T 



2 2 2 

m a ~ m ~ m 7T 



aor)NTr m a J 



!7r (m) = A a 



(B5) 
(B6) 



and 



A aoVN TT = Zl(d(j) 2 s - A 2 )^JV, 



B, 



a i] N TT 



)9 2 ^ N \^ IZlf 

2 TO 2 



N 



{hi - h 3 ) 



(B7) 
(B8) 
(B9) 



4 Obviously, the neutral and the charged ao's have the same formulas for the decay widths. 
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AaorisTr — -^ClZTrZ ris <t> N <t)S- 



(BIO) 



The ao — > KK decay width includes two subchannels, K°K° and K K + , which results in the isospin factor 1=2, 
and accordingly the decay width is found to be 



Tan 



* /i-f^Y 

8ro ao 7r y \ m ao J 



A ao KK - -^B aa KK{m 2 aa - 2m 2 K ) + C ao KKm 



where 



Aa 



KK 



\ 2 Z\ \ 9N 



4>S ' 

V2j' 



B aa KK = Z 2 K w Kl j.gi - -wk x {{gl + h 2 )(f>N + V2(gj - /i 3 )</>s) j , 



(BH) 

(B12) 

(B13) 
(B14) 



Now turning to the scalar kaon there is only one relevant decay channel, Kq (or K ) — > Kir, for which the decay width 
reads 



Kq — yK-ii 



8irm Ko 



("4-q - m l - m k) 2 - 4774774-1 

4ml 

Kn 



Ak Kit + {CkoK-k + DkqKtt — Bk Ktt) 



l K 



CkqK-k^k + B>K KTrm v 



(B15) 



where 



Ak Ktt 
BkqKtt 



-A 2 <Ps, 



V2 
Z-kZkZk 



-W ai WK 1 



2gx 



w ai + W Kl 
w ai w Kl 



+ (2h 3 -h 2 - Zg{)4>N - Vz(gl + h 2 )<j> s ) 



CkoKtt = * K ° [~gi{iwK* + w Kl ) + V2iw K *w Kl {g 2 ~ h 3)^s], 

D Ko Ktx = — ^— ; — — [1gi{iwK* ~w ai ) + iw K *w ai ((2h a - h 2 - 3gj)(f> N + V2(gj + h 2 )<j>s)- 



(B16) 

(B17) 

(B18) 
(B19) 



Finally, for the f particles there are two relevant decay channels: the two-pion and the two-kaon channels. It is 
important to note that due to the particle mixing between and f$ , the matrix elements are linear combinations 
that depend on the scalar mixing angle a (see Appendix [3} , as can be seen explicitly below. The decay widths in 
the irir channel are 



J ° 327TTOyL \ 



2m, 



771 x L 
* ^0 , 



M 



r fH _ 



2m, 



\ \ m / H 

where an isospin factor of 3/2 was considered 5 and the matrix elements are 
M f t = - sin V M% V (m f L ) + cos O a M^( m tf)> 



(B20) 
(B21) 

(B22) 



5 There are two subchannels, namely the n + ir~ and the 7r°7T , which would mean I = 2. However, since the two 7r° are indistinguishable, 
the solid-angle integral (there is no angular dependence at tree-level) ends up as 2w instead of Air, which means a factor of 1/2 in case 
of 7r 7r° compared to 7r+7r _ , thus I = 1/2 + 1. 
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M^(m) = 2Zl4>N 
Mf on (m)=2Z*ct>si 



g| m 2 r / _ 2m^ ?7i 2 + /ii0|/2 + 2(5^ 

2 < 

gi "z 2 
4 m 2 ai 



1 - 



2m 2 \ hx(j) 2 N 



-Ai 



In the if If channel, where 1 = 2, the decay widths read 

1 



r /«^AA 



fcf-tKK 



1 



2m,K 



1 



87rm 



2m,K 
irifL 
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where the matrix elements, using the notations i/jv = (<?? + 2/t.i + /12) /4 and Hs = (g\ + h\ + /12) /2, are 



M L faK {m) 



Mf oK (m) 



- sin e a Mf oK (m f L ) + cos 6 a M% K (m ^ ) 
cos 0aM% K (m f H ) + sin 6> CT M^ ^ (m /( h ) , 
A 2 



Z A 



(2A 1 +A 2 )0 A r 

2(Ai + A 2 )0 s - 



^PS + 9iwk 1 [m K - m ) + w Kl \ 2H N (/) N ^f^^ 5 
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A2 
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<p N + V2g 1 w Kl (m 2 K - m 2 ) 



w 



I 2Hr^ h3 ~ 9 ^-^ m2 ~ 2m ' K 



K 



ISPS 



V2 
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2. Vector-meson decay widths 

In the case of the vector mesons we are considering the decays of the p meson, the K* vector kaon, and the <f> meson. 
All these particles have only one relevant decay channel. The first is the p — !> nn decay which has the following quite 
simple decay width formula 



p—>7t7T 



487T7< 



/ 2m w 



3/2 



rj2 92 , rj2 



The next one is the K* — > Ktt decay, in case of which the decay width reads 



A*-s-Aj 



m K * 
8tt 



A" 



4m|, 



3/2 



(A k * K tt — B K * Kyr + Ck*k-it) 2 



where the constants are defined as 
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A K* Ktt = -Z^Zr 
1 

Bk'Ktt — -^Z^Zk 



-gi + v2w Kl (gf - h 3 )<f>s 
2.9i + w ai (-3gj -h 2 + 2h 3 )(p N + V~2w ai (g\ + h 2 )(j>s 



Ck*Ktt = ■^Z^Z K w ai w Kl g2 , m 1 K* 



Finally, the <f> — > KK decay width reads 



192nm 4 Ki 



1 - 



/ 2mA" 
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-I 2 
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3. Axial-vector-meson decay widths 



Turning to axial-vector mesons, the considered decays are the a\ decay with two relevant channels and the /is 
decay with one relevant channel. Since in case of the a\ — > pit and fig — > KK* decays, the decaying as well as one 
of the resulting particle are (axial-)vector mesons, the matrix elements have a more complicated form than in the 
previous cases, as can be seen below. The first considered decay width is the one of the ai — ► TT-f process, which takes 
the following simple form 



e 2 gl4> 2 N Z 2 

967rm ai n 



m„. 



(B36) 



In case of the a\ — > pn decay width one has to consider two channels (p + n , p tt + ), thus 1 = 2 and the decay width 
is found to be 



1 



ai —>p-K 



I2m ai ir 



(m„j ~ m 2 — m 2 ) 2 — 4m 2 m 2 



1/2 



m 2 p m 2 ai 
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where Vu V is the vertex following from the relevant part of the Lagrangian, 



= iZ^(j) N I (g\ - h 3 ) + -^-[k^^ka^ + k pp k nu - k n ■ (k p + k ai )g pv 



(B38) 



and fc^ = (m ai ,0), kp 1 = (E p ,\c) and k% = (E n ,—li) are the four- momenta of oi, p, and tt in the rest frame of a\ 1 
respectively. Using the following kinematic relations 



k ■ h 



™> 2 ai - m l- m l 
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k ai ' kp — m ai E p 
k 2 = 



m ai + mf p -m$ 



(m 2 ai -ml~ m 2 ) 2 - ^m\m 2 



the terms in Eq. (|B37|) are given by 
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Analogously to the previous case, the width of the fis — > KK* decay which includes four sub-channels (A"°A"*°, 
K°K*°, K~K*+, K+K*~), resulting in 2 = 4, becomes 
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>KK* 
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4m 4 f 
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where V pu has the same Lorentz structure as Eq. (|B38|) and the only difference is in the constants in front of the 
different terms. More explicitly, 



Vpv = iZ K {A fKK g pv + B fKK [k Kp kf lsl , + k K * p k Ku - k K ■ (k K * + kf is )g pu ]} , 



(B42) 
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A f KK = t \gl(V24> N - 6(/>s) + V2h 2 (<j) N - y/2</> s ) + 4/i 3 & 



B f KK = --j=g 2 w Kl - 



(B43) 
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The kinematic relations are the same as in Eq. (jB39| with the following substitutions: ai — > fis,ir — > K, p — > K*, 
while the expressions analogous to Eq. (|B40[) are, 



\V P „\ 2 = ZH4A 2 firir + B 2 f 



K I }KK ~ ^fKK 
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I V^kfr I 2 = ^ {A 2 kk B 2 KK {E 2 K * -m 2 K *) + 2A fKK B fKK (E 2 K - m 2 K )} , 
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2 2 

K lis 



A 2 

7 1 fKK ^2 
6 K 2 ^K* ■ 

m K* 



Analogously, the width for a generic decay of the form K\ — > VP, where V denotes a vector and P a pseudoscalar 
state, reads 



K^VP 



= 1 



k(m,K 11 m,v ,mp) 1 
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where I = 3 for the Ki —> pK and A'! — > A'*7r decays and X = 1 for the A x — > wjvA" decay, 

fc(m , 1716, ™ c ) = tt" - J - 2ml ( m b + m c) + ( m l ~ miy 2 8(m a -m b ~ m c ) 
2m a V 

(the theta function ensures that the decay width vanishes below threshold) and 

{A Kl g^ + B Kl [p VPp + PpP » Ki - ( PP ■ p v )g^ - (p Kl ■ p P )g» v ] } . 



v KiVP 



(B47) 
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The values of Ak %1 Bk x1 and Ck x depend on the process considered: Ak x — {Ak 1 k*tt, Ax iP k, Ak iU!n k} and Bk 1 
{B KlK * 7r ,B KlpKl B KlUlNK }. The coefficients read 



A Ki k*ti = —j=Z^(h z - g\)(j>s , 
BkiK*7c = —-^Z 7l g 2 w ai 



Ak iP k = 
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Ak iUn k = -^Z K [g\{4> N + V2(f>s) - h 2 (<p N - V2</>s) - 2h 3 <f) N ] , 

BKna N K = —-jZKg2WK x ■ 
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